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I  EXECUTIVE  SUMMARY 


This  report  completes  the  propagation  modeling  work  performed  under 
DNA  Contract  DNA001-77-C-0038  (SRI  Project  5960).  The  broad  objective 
of  this  work  was  to  develop  a  unifying  theoretical  framework  for  the 
various  nuclear- propagation  predictions  based  on  (1)  numerical  simulations, 
(2)  field  measurements,  using  the  naturally  and  artificially  disturbed 
ionosphere,  and  (3)  predictive  computer  codes  for  systems  analysis. 

It  is  now  generally  accepted  that  knowledge  o:  the  striation 
spectral-density  function  (SDF)  suffices  for  pret  -ing  the  signal 
moments  that  characterize  the  average  disturbed-.. ignal  structure.  The 
mathematical  form  of  the  SDF  follows  a  power  law.  The  range  of  the  power- 
law  continuum  and  the  appropriate  spectral  index,  however,  admit  some 
uncertainties . 

Naturally  occurring  phase  scintillation  admits  a  power- law  spectral 
characterization  over  a  scale-size  range  from  tens  of  kilometers  to 
hundreds  of  meters  with  a  one- dimensional  power- lav  index  near  2.5.  For 
typical  barium  releases,  the  power- law  continuum  evi’intly  does  not  extend 
to  scale  sizes  larger  than  one  kilometer.  The  corresponding  phase  spectrum, 
derived  from  back- propagated  LES-8  signals  during  the  STRESS  experiments, 
gives  a  power- law  index  somewhat  greater  than  three  [Knepp  and  Bogusch, 
1979].  These  results  are  not  necessarily  in  conflict,  but  they  do  imply 
that  the  propagation  theory  must  accommodate  some  variation  in  spectral 
slope. 

The  theoretical  work  performed  under  this  contract  fully  accommodates 
a  varying  spectral  index.  The  results  show  that  the  signal  structure  in 
a  power- law  environment  is  critically  dependent  on  the  power- law  slope, 
particularly  under  strong-scatter  conditions.  The  effects  are  subtle, 
but  they  do  affect  predictive  modeling. 
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The  predictive  modeling  problem  is  typically  broken  into  two  parts. 

The  weak-scatter  theory  is  linearly  extrapolated  to  determine  the  onset 
of  saturated  scintillation.  The  Rayleigh  model,  which  is  completely 
characterized  by  the  mutual  coherence  function,  is  then  assumed  to  apply. 
This  scheme  is  generally  adequate.  Our  work  shows,  however,  that  in  a 
power-law  environment  the  Rayleigh  limiting  form  can  only  be  strictly 
achieved  when  the  phase  spectral  index  is  less  than  3,  as  is  evidently 
the  case  for  naturally  occurring  scintillation. 

It  is  intended  that  this  report  be  self  contained;  thus,  some  back¬ 
ground  material  is  presented  in  Section  II.  In  Section  III  we  summarize 
the  we.ak-scatter  theory  that  characterizes  the  onset  of  amplitude  and 
phase  scintillation.  A  detailed  development  can  be  found  in  Rino  (1979a). 
In  Section  IV  the  general  formulas  for  the  phase-screen  model  are 
developed.  The  general  formulation  fully  allows  for  wavefront  curvature, 
which  was  not  explicitly  included  in  our  previous  calculations  based  on 
incident  plane  waves.  The  spherical-wave  correction  factors  are  easily 
applied  to  the  previously  developed  plane-wave  results. 

In  Section  V  a  model  for  the  mutual- coherence  function  is  developed 
and  then  used  to  derive  simple  formulas  for  the  signal- coherence  times. 

Two  approximations  are  used:  one  valid  for  shallowly  sloped  spectra;  the 
other  for  the  more  steeply  sloped  spectra.  In  Section  VI  a  frequency- 
coherence  model  is  developed  and  simple  formulas  for  the  random  delay 
jitter  and  pulse  dispersion  derived  from  it. 

Finally,  in  Section  VII  a  model  for  angle  scatter  is  developed.  In 
all  cases  the  formulas  fully  accommodate  propagation  at  arbitrary  incidence 
angles  relative  to  the  magnetic  field.  In  deriving  several  of  the  for¬ 
mulas  that  might  ultimately  be  used  in  a  predictive  code,  we  used  a 
quadratic  approximation,  which  is  not  strictly  valid  for  the  shallowly 
sloped  spectra.  Some  refinements,  therefore,  may  ultimately  be  desirable. 

Wideband  satellite  data  as  well  as  data  from  the  1979  PLUMEX  rocket 
campaign  at  Kwajalein  are  currently  being  analyzed  and  '-.hould  show  the 
degree  to  which  the  approximations  are  valid,  at  least  in  the  natural 
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environment.  In  their  present  form,  or  with  rather  small  modifications, 
the  formulas  sunmarized  herein  are  adequate  for  current  predictive 
modeling  requirements. 


II  INTRODUCTION  AND  BACKGROUND 


One  objective  of  our  channel  modeling  work  has  been  to  develop  a 
mathematically  tractable,  but  accurate  characterization  of  the  onset  of 
phase  and  amplitude  scintillation  and  the  temporal,  spatial,  and  frequency 
coherence  of  the  disturbed  radiowave.  The  coherence  measures  can  be 
directly  related  to  the  scincillation-induced  degradation  of  a  variety 
of  systems. 

We  have  pursued  a  formal  channel  modeling  approach,  which  takes 
advantage  of  the  linearity  of  the  transionospheric  channel.  If  the 
complex  signal,  v(t)s  is  transmitted,  fpr  example,  the  received  signal, 
vQ(t),  [less  noise  and  interference!  admits  the  representation 

vQ(t)  =  J v(f)h(t;f+fc)  exp  {2rtift}  df  (2.1; 

A 

where  v(f)  is  the  Fourier  transform  of  v(t),  f£  is  the  center  frequency, 
and  h(t;f)  is  the  time- varying  transfer  function  that  characterizes  the 
channel . 

Formally,  h(t;fc)  is  the  response  of  the  channel  to  a  sinusoidal 
signal  v(f)  ~  §(f~fc).  The  time  structure  comes  about  primarily  because 
of  the  relative  motion  of  the  propagation  path  within  the  randomly 
irregular  medium.  In  effect,  the  receiver  sees  a  moving,  spatially 
modulated,  spherical  wavefront.  If  u(p,z;fc)  represents  the  spatial 
modulation,  then 

h(t,fc)  =  u(vTt,vzt;fc)  (2.2) 

where  v  =  (v  ,v  )  is  the  instantaneous  relative  velocity.  The  3cintil- 
lation  problem,  therefore,  is  mainly  one  of  characterizing  the  spatial 
modulation  imparted  to  a  spherical  wave  by  a  random  medium.  The  time 
structure  is  obtained  by  a  straightforward  translation. 
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There  are  two  time  and  associated  frequency  intervals  of  interest 

in  all  communication  systems:  (1)  a  comparatively  short  time  interval 

and  its  associated  instantaneous  bandwidth  about  f  ,  which  determines  the 

c 

basic  modulation  scheme  (e.g.,  FSK,  DPSK) ;  and  (2)  a  much  longer  time  and 

frequency  interval,  over  which  coding  and  frequency  hopping  are  employed. 

In  almost  all  systems,  h(t;f+f  )  is  nearly  invariant  over  the  sub- frequency 

c 

band  occupied  by  v(f).  It  them  follows  from  Eq.  (2.1)  that 

v  (t)  —  h(t;f  )v(t)  .  (2.3) 

O  v 

Scintillation  is  simply  a  complex  (amplitude  and  phase)  modulation 
imparted  to  the  signal.  Going  further,  if  the  time  scale  of  h(t;fc)  is 
long  compared  to  the  duration  of  v(t),  then  h(t;f£)  is  effectively  a 
complex  constant.  In  such  a  "slow-fading"  environment,  only  the  amplitude 
variations  affect  the  performance  of  a  well-designed  system. 

As  a  practical  matter,  the  Nakagami  distribution 

m  m- 1 

P(y)  =  mr(m~)~  exP  f-my}  (2.4) 

_2 

where  y  is  the  signal  intensity,  m  =  is  the  intensity  scintil¬ 

lation  index,  and 

S l  -  «I2>  -  <l)2)l/2/<I>  ,  (2.3) 

gives  a  very  good  approximation  to  the  intensity  statistics  under  all 
scatter  conditions  [see  Fremouw  et  al.  (1980)].  Moreover,  under  strong- 
scatter  conditions,  ~  1  so  that  m~  1  and  Eq.  (2.4)  takes  the  exponen¬ 
tial  form  of  the  Rayleigh  distribution. 

We  now  know,  however,  that  the  complete  Rayleigh  model  is  not 
strictly  applicable  to  signals  undergoing  strong  scattering  because  of 
omnipresent,  slow  phase  and  associated  amplitude  variations.  The  impact 
of  such  subtle  departures  from  strict  Rayleigh  scattering  is  discussed  in 


Rino  (1979a, b) .  For  predictive  modeling,  the  departures  from  predictions 
based  on  the  strict  Rayleigh  model  are  generally  negligible. 

The  power  of  the  Rayleigh  model  is  that  it  is  completely  specified 
by  the  mutual  coherence  function 

Ru(8t)  -  <h(t;fc)h*(t';fc»  ,  (2.6) 

which  is  easily  computed.  In  fact,  for  many  applications  one  need  only 
specify  the  coherence  time,  t  ,  such  that 

VTc>  -  K  V0)  (2-7> 

where  K  «  e  ^  or  some  other  convenient  fraction.  In  all  our  applicati^w?. 
R^O)  ■  (i)  ■  1,  that  is  the  signal  is  normalized  to  unit  intensity.  A 
second  useful  relationship  for  Rayleigh  channels  is 

<l(t)I(t'»  -  1  -  |Ru(6t)|2  .  (2.8) 

It  follows  from  Eq„  (2.8)  that  the  intensity  coherence  time,  Tj,  and  the 

coherence  time  of  the  complex  signal,  t  ,  are  simply  related. 

c 

The  first  step  in  predictive  modeling,  therefore,  is  to  determine 
the  onset  of  saturated  scintillation  (S^  ~  1),  where  the  Rayleigh  model 
can  be  applied.  The  system  effects  under  weak  scintillation  are  easily 
determined  from  Eq.  (2.4)  and  a  rough  estimate  of  the  fade  coherence 
time.  To  estimate  the  onset  of  saturated  scintillation,  wc  linearly 
extrapolate  the  weak  scatter  theory  beyond  its  strict  range  of  validity. 
This  is  a  very  conservative  estimate  because  strong  scatter  effects 
always  act  to  reduce  the  weak-scatter  value.  A  detailed  treatment 
of  strong-scatter  effects  is  presented  in  Rino  (1979b,  1980).  The  weak- 
scatter  theory  is  developed  in  Rino  (1979a),  but  for  convenience  the 
principal  formulas  are  summarized  in  Section  III. 
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All  our  theoretical  developments  are  based  on  the  phase-screer 
model  not  only  because  of  its  simplicity,  but  because  recent  ana'  .s 
show  that  results  based  on  the  phase-screen  model  preserve  all  essential 
aspects  of  much  more  elaborate  calculations  that  properly  accommodate 
the  effects  of  diffraction  within  .'he  scattering  medium  [Bramley  (1977)  ; 
Fante  (1976)].  Furthermore,  results  based  on  the  phase- screen  model 
can,  in  many  cases,  be  applied  incrementally  to  accommodate  variation 
along  the  propagation  path  (Rino,  1978). 

In  all  our  previous  analysis,  wa  have  treated  only  the  scattering 
of  an  incident  plane  wave.  The  plane-wave  results,  however,'  are  easily 
generalized  to  accommodate  the  incident  spherical  wavefronts  that 
eminate  from  a  source  at  a  finite  distance.  This  is  demonstrated  in 
Section  IV  where  the  phase-screen  model  is  applied  to  calculate  the 
mutual  coherence  function. 

The  multiple-scatter  theory,  based  on  solutions  to  the  parabolic 
wave  equation,  depends  only  on  the  integrated- phase  structure  function, 
as  does  the  phase-screen  theory.  Because  of  its  central  role  in  the 
propagation  theory,  the  form  of  the  phase  structure  function  in  a  general 
pc;wer- law  environment  and  various  approximations  to  it  were  discussed 
In  detail  In  Topical  Report  4.  In  Section  V  of  this  report  we  use  those 
results  to  compute  the  mutual  coherence  function  Eq.  (2.6). 

Precision  navigation  systems  use  very  large  bandwidths,  and  questions 
of  frequency  coherence  become  important.  Thus,  we  have  also  computed 
the  appropriate  frequency- coherence  function  and  used  it  to  estimate 
the  average  delay  jitter  and  coherence-bandwidth  loss.  The  calculations 
are  summarized  in  Section  VI.  To  complete  the  channel  model,  in  Section 
VII  we  compute  the  angular  spectrum  and  derive  formulas  that  characterize 
the  angle  scintillation.  These  formulas  can  be  used  for  radar  appli¬ 
cations  as  well  as  relating  the  model  calculations  to  previously  developed 
radar  codes. 

In  all  cases,  the  model  calculations  reflect  the  critical  dependence 
of  the  scintillation  structure  on  the  power" law  index.  If  the  index  of 
the  one-dimensional  phase  spectral  density  function  is  less  than  3,  the 
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statistical  theory  can  be  formulated  completely  without  specifying  the 
outer  scale.  As  the  spectrum  steepens  so  that  the  phase  spectral  index 
exceeds  3,  however,  the  outer  scale  must  be  specified  to  compute  the 
mutual  coherence  function.  In  our  model  development  we  have  generated 
formulas  that  apply  to  Loth  the  shallowly  sloped  and  steeply  sloped 
spectra. 
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Ill  THE  WEAK- SCATTER  THEORY 


For  modeling  purposes,  we  can  ignore  the  inner  scale  cutoff.  This 
has  no  effect  on  the  phase  variance,  and  at  worst  it  over  estimates  the 
intensity  scintillation  at  the  highest  frequencies  of  interest.  The 
three-dimensional  SDF  that  characterizes  the  striations  has  the  form 


abC 

$£N  ^  =  2.V+1/2 

e  [q  +  q  1 


(3.1) 


where  q  is  the  outer- scale  wave  number,  and  a  and  b  are  the  axial  ratios 
o 

along  and  transverse  to  the  magnetic  field  respectively.  The  morr'. 
general  form  of  Eq.  (3.1)  that  includes  an  inner-scale  cutoff  is  discussed 
in  Topical  Report  4. 

The  isotropic  turbulent  strength  C  is  related  to  the  rms  electron 

2  s 

density  (ANe)  by  the  relation 

Cs  =  8rt3/2  <M2)  q2v"2  r(vfl/2)/r(v-l)  .  (3.2) 

2 

Since  (£N  )  and  q  or  the  outer  scale  1  ~  2«/q  are  most  often  specified, 

CO  0  0 

Eq.  k'3.2)  can  be  used  to  calculate  Cg.  The  autocorrelation  function  of 
the  integrated  phase  has  the  fomi 


(y)  =  r  2x24  G  C 
e  p  s 


v‘1/2  Vi/2(y> 

2rtT(vfl/2) 


(3.3) 


where  rg  is  the  classical  electron  radius,  \  is  the  wavelength,  is 
the  length  of  the  propagation  path,  and  ^(x)  is  the  modified  Bessel 
function. 

To  account  fully  for  the  spatial  variation  of  the  integrated  phase, 
y  in  Eq.  (3.3)  must  be  replaced  by 


£</*.*)  - 


2  211/2 

C*P8  -  BZ^s  Apg  +  AAp^  U 

x _ x  y _ y 

AC  -  B2/4 


(3.4) 


where 


4pc  =  As*  -  tan  0  a.  Az 


(3.5) 


The  angle,  0,  is  the  angle  of  the  propagation  vector  relative  to  the  z 

A 

axis  and  a^  is  a  unit  vector  along  the  projection  of  the  propagation 
vector  onto  the  plane  normal  to  z.  The  reason  Eq.  (3.5)  takes  this 
particular  form  will  become  clear  in  Section  IV. 

The  coefficients  A,  B,  and  C  depend  on  the  relative  propagation 
angles  and  the  anisotropy  model.  The  geometrical  enhancement  factor,  G, 
is  defined  in  terras  of  A,  B,  and  C  as 


(3.6) 


/AC  -  B  /4 


The  anisotropy  model  is  summarized  in  the  Appendix. 

The  temporal  phase  autocorrelation  function  is  obtained  by  sub¬ 
stituting  v,p6t  and  v  6t  for  Zp  and  Az  respectively  in  Eq.  (3.5).  The 
1  2  X 

final  form  can  be  written  compactly  as  where 


Cv2  -  Bv  v  +  Av2  11/2 
s  s  s  s 

X  XV  v 


AC  -  B  /4 


(3.8) 


vg  =  vT  -  tan  0  a^vz 


(3.8) 


The  velocity  vector  v  =  (v„,v  )  is  formally  the  relative  velocity  of  the 

1  z 

irregularities  as  seen  by  an  observer  in  the  reference  coordinate  system. 
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The  one- dimensional  phase  SDF 
transform  (veff^  *■)  *  The  result 


is  obtained  by  taking  the  Fourier 
is 


(p(£) 


[f2  +  f2]p/2 


(3.9) 


where 

fo  =  veff  qo/2rt  =  veff  lo  >  (3*10) 

and 


T 


IGC  - - 

P  S  (2rt)2v+V(v+l/2) 


2v-l 


v 


eff 


(3.11) 


In  Eq.  (3.9),  p  =  2v. 

The  phase  variance  can  be  computed  by  integrating  Eq.  (3.9)  over 
all  temporal  frequencies,  f.  The  result  is 


2  2  2  qo2V+1  r<v-1/2> 

(60  >  =  reX  *pG  Cg  4rt  r(v+l/2)  * 


In  most  cases,  however,  the  length  of  the  data  interval,  t  ,  is  less 

-1  c 
than  f  so  that 
o 


<602) 


2T 

p-1 


P-1 


(3.13) 


To  estimate  the  measured  rms  phase,  one  should  use  the  minimum  value  of 
Eq.  (3.12)  or  Eq.  (3.13). 


The  corresponding  intensity  s cintillation  is  given  by  the  formula 
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(3.14) 

(3.15) 

(3.16a) 

(3.16b) 

(3.16c) 

(3.17) 


where  is  the  signal- source- to- phase- screen  distance  and  is  the 
phase-screen-to-receiver  distance.  A  simple  series  expansion  can  be 
used  to  evaluate  the  hypergeometric  function.  For  actual  computations 
we  take  the  minimum  value  of  Eq.  (3.14)  or  unity.  A  predicted  value  of 
unity  from  Eq.  (3.14)  is  a  conservative  estimate  of  where  saturation 
actually  occurs. 

Equation  (3.14)  is  strictly  valid  only  when  \/ZqQ  «  1,  which  is  a 

good  approximation  at  or  above  VHF  for  most  propagation  conditions.  If 

*/Zqo  »  1,  however,  it  is  easily  shown  that  ~  2ct^.  If  «s/Zqo  >  1, 

the  limiting  value  rf  2o^  can  be  used  in  place  of  the  value  given  by 

Eq.  (3.4)  as  a  simple  means  of  accommodating  far-zone  scattering.  In 

an  extended  medium  l  can  be  treated  as  a  differential  distance  and 

P 
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Eq.  (3.12)  or  Eq.  (3.13)  and  Eq.  (3.14)  integrated  along  the  propagation 
path  to  accommodate  the  varying  structure  and  propagation  geometry. 

Once  the  signal  has  saturated,  one  ca.'  use  the  Rayleigh  model  to 
predict  the  signal  structure.  As  discussed  in  Section  II,  however,  the 
strict  applicability  of  the  Rayleigh  model  is  open  to  question.  The 
power  of  the  model  is  that  it  is  completely  specified  by  the  mutual 
coherence  function.  We  also  note  in  passing  that  under  weak-scatter 
conditions  the  intensity  coherence  time  generally  satisfies  the 
inequality, 

TI<v'z/veff  .  (3.18) 


This  is  discussed  in  detail  by  Rino  and  Owen  (1980) 


IV  THE  PHASE- SCREEN  MODEL 

The  phase-screen  model  proceeds  formally  from  the  Huygens-Fresnel 
principle,  which  for  our  purposes  can  be  written  as 


u(r)  =  C 


H 


u(r')G(r,r')  dr' 


where 


exp  f-ik  r  -  r'  } 


r  -  r 


(4.1) 


(4.2) 


and  the  complex  constant  C„  will  be  determined  later  [see  Born  and  Wolf 

n 

(1964),  Chapter  8  for  a  detailed  discussion].  If  u(p,zQ)  is  known,  we  can 
write 


u(p,z)  =  C 


H 


ff 


rr 


Oc;zo) 


G(r,r/)  exp  f-i^  •  (p  -  p ')}  dp' 


X  exp  f-iK  •  p} 


die 

(2rf)2 


(4.3) 


where  u(K;zq)  is  the  Fourier  spectrum  of  u(p,ZQ).  For  notational  con¬ 
venience  we  let  zq  =  0  when  its  value  enters  explicitly. 

The  inner  integral  can  be  evaluated  [see  Rino  and  Fremouw  (1977), 

Eq.  (5)].  Substituting  that  result  into  Eq.  (4.3)  and  changing  variables 

— ♦  *4 

from  ic  to  K  +  k^,  giyes 


4' 


u(p,z)  =  -iXC 

X  exp  f-ik^,  •  p] 


_  .  exp  f-ikg((?  +  lL)z}  _ 

(K  +  kj)  - - - - -  exp  {-iK 


Pi 


dtf 


g(K  +  Kj.) 


(2*)‘ 


(4.4) 
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We  have  set  C  =  i  sec  q/\  so  that  Eq.  (4,8)  is  a  proper  solution  to  the 
parabolic  wave  equation. 

Now  consider  an  incident  spherical  wave 
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ui = 


exp  J  -ik|Rx  +  p|j 
Rx  +  p[ 


(4.11) 


where  R.^  is  the  distance  from  the  point  source  to  the  plane  z 
p  lies  within  the  plane.  If  »  p,  we  can  write 


z  and 
o 


(4.12) 


Replacing  the  reference  wave,  exp  {ikR^}/Rp  by  a  more  general  modulation, 


u<p,z  ),  we  have 


u(K  +  k) 


■J}"®' 


2  2 

p7  -  sin  0 


,z  )  exp  <- 


(V?< 


Ri 


Substituting  Eq.  (4.13)  into  Eq.  (4.8)  gives 


exp  {-iK  •  p7}  dp7 
(4.13) 


2 


u(p,z) 


p7  -  sin  0 


u(p7,zq)  exp  ^ 


(\  ' p')' 


R. 


i  V  )  ' 

exp  |-ih(K)  -  1  exp  {iK  •  (p7  -  pg)} 


A-  &  (4.14) 

(2rt) 


The  integral  over  K  can  be  evaluated  as 


T  j  cos  0  /  j  m 

1  =  -1  exp  " 


n2  -  sin2  0(akT  .  Ti)2) 


z  sec  0 


i 


where  T]  =  p  -  pg  and  R£  is  the  propagation  distance  from  z 
(R2  =  z  sec  0) .  Thus 


=  z 


(4.15) 


Y 
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.  lc  cos  0 

u(p’z)  =-2^-])u(p> 


Zo}  eXp  p  2 F 


(p#  -  Cp, 


v  2  ,  2«  /a  — * / 

)  -  sm  P  " 


Cpr 


df 


X  exp{-i 


<R1  +  V 


p~  -  sin“e^  •  P„\ 


>2S  --2' 


-7 


(4.16) 


where 


R_  =  c: 


R1R2 


e  Rx+R2 


(4.17) 


and 


ti. 


C  = 


Rl  +  R2 


(4.18) 


It  follows  from  Eq.  (4.16)  that  aside  from  a  small  phase  modulation,  the. 
effects  of  a  point  source  of  illumination  can  be  accommodated  simply  by 
replacing  pg  by  Cpg  and  the  propagation  distance  (z  sec  8)  by  Rj^/ (R^+R,,) 
Thus,  any  result  based  on  plane  waves  is  readily  modified  to  accommodate 
spherical  waves. 

Equation  (4.16)  also  shows  the  effects  of  an  interchange  of  source 
and  receiver  (R^  and  R2) •  Insofar  as  effects  that  depend  on  propagation 
distance  are  concerned,  reciprocity  holds.  Structure  variations,  however, 
depend  on  C,  and  they  do  change.  As  a  practical  matter,  only  and  other 
single- point  measurements,  such  as  the  frequency- coherence  measures 
described  in  Section  VI  are  reciprocal.  Spatial  and  temporal  coherence 
measures  (second-order  moments)  are  not  reciprocal. 

Equation  (4.16)  or  its  Fourier  domain  equivalent  form  the  starting 
point  for  all  phase-screen  computations.  Hereafter,  however,  we  shall 
apply  the  plane  wave  form  and  simply  insert  the  spherical  wave  correction 
factors  where  appropriate.  This  has  already  been  done  in  Eq.  (3.17) 


V  THE  MUTUAL  COHERENCE  FUNCTION 


/ 


As  discussed  in  Section  II,  under  conditions  of  Rayleigh  scattering 
the  mutual  coherence  function  completely  specifies  the  average  sign'll 
structure.  To  compute  (u(p,z)u*(p;z ')) ,  we  use  the  Huygens-Fresnel 
formula  Eq.  (4.8)  together  with  the  homogeneity  assumption 

<u(j?)u*(K'))  =  $u(K) S(K  -  K')  (5.1) 

as  discussed  by  Rino  and  Fremouw  (1977).  In  Eq.  (5.1),  $U(K)  is  the 
SDF  of  u(p,ZQ)  and  6(K)  denotes  the  Kronecker  delta  function.  The 
result  is 


<u(p,z)u*(p'z')>  =  |  $u("ic  +  Kj)  expj-ih(K) 


Az  sec  Q i 


X  exp  J-iK  •  Ap 
( 


die 

(2  7i)2 


(5.2) 


where  Az  =  z  -  z',  Ap  =  p  -  p',  and  Apg  -  Ap  -  tan  0  a^Az, 

It  is  simpler  to  specify  the  mutual  coherence  function  than  the 
angular  spectral  density  §U(K  +  kj) .  We  note,  however,  that  if  Az  =  0, 
Eq.  (5.2)  shows  that 


<u(p,z)u*(p',z)>  =  RuC*ps) 


(5.3) 


Indeed,  for  most  applications  the  Az  dependence  can  be  ignored.  Nonethe¬ 
less,  by  using  Eq.  (4.15)  it  is  readily  shown  that 


Ru(4p,Az)  = 


jj 


„  S  k  cos  9 

V *  +  *>•>  2^~ 


X  exp<-i 


2^  [^'2  ’  SI"2  6(\  '  ^')2]|<^'  .  (5.4) 


where  we  have  substituted  for  Az  sec  0. 
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In  the  phase-screen  models 


Ru0£)  =  exp  {-  |d^(^s)} 


(5.5) 


where  Dg^(Zip)  is  the  integrated- phase  structure  function.  Two  approximate 

forms  for  D  (y)  have  been  developed.  If  the  v  parameter  [see  Eq.  (3.1)] 

o(»  " 

is  less  than  1.5, 

(5.6) 


Vy) 


c2  |y|2v“1 
60m 


where  C  is  the  phase- structure  constant 


C2  =  G  — ^ 
C6 <t>  G  2rt 


zrd.  5-v) 


.  2v-l 


(5.7) 


2,2 


In  Eq.  (5.7),  C  .  rVj l  C  , 
’  p  e  p  s 


r(v+0.5)(2v-l)‘ 

A  detailed  derivation  of  Eq.  (5.6)  and  a 


demonstration  of  its  range  of  validity  can  be  found  in  Topical  Report  4. 

If  v  2  1.5,  the  situation  is  more  complicated.  The  overall  conver¬ 
gence  of  the  Taylor  series  is  too  slow  for  it  to  be  of  general  practical 
applicability.  It  happens,  however,  that  the  quadratic  term  alone  does 
an  acceptable  job  of  reproducing  the  exact  form  of  Dg^(y)  f°r  %|y|  «  1 
when  v  >  1.5.  Thus,  we  take 


Vy>  ~  *v”i<v>2 


(5.8) 


where 


D 


1 


!\  log  (2e) 

r(v-1.5) 

4r(v-0.5) 


v  =  1.5 


v  >  1.5 


(5.9) 


and  e  «  1  is  the  ratio  of  the  outer-  to  inner-scale  wave  numbers.  Note 
that  Eq.  (5.8)  depends  on  qQ  and,  when  v  =  1.5,  the  inner  scale  as  well. 
Again,  a  detailed  ..iscussion  is  given  in  Topical  Report  4. 
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In  general,  therefore. 


%<y)  ~  Dv|y|min(2v"1»2)  (5.10) 

2  2 

where  D  is  replaced  by  Eq.  (5.7)  when  v  <  1.5  and  2o\D,q  with  D,  given 

910  1 

by  Eq.  (5.9)  when  v  £  1.5.  Unfortunately,  Eq.  (5.4)  cannot  be  evaluated 
analytically  when  v  <  1.5.  If  v  s  1.5,  however,  it  can  be  shown  that 


(5.11) 

If  |d£>|  »  AR2  as  is  usually  the  case,  then 

R^,^)  ~  expj-  |  Dv|^s|min(2v"1’2>J  .  (5.12) 

If  »  | A? I,  however, 

|r u(z£,ar2)  I  ~ 


^2? 
k  / 


provided  that  8  is  not  too  large. 

To  characterize  the  measured  signal  when  Eq.  (5.12)  applies,  we 

have 

V'O  =  expj-  \  DvJcveff6t]min(2v''1»2)|  .  (5.14) 
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The  C  factor  as  defined  by  Eq.  (4.18)  corrects  for  wavefront  curvature. 
If  Eq.  (5.3)  applies,  we  have,  similarly, 


/e  hn 


where  vm  is  the  rate  of  change  of  R  as  defined  by  Eq.  (4.17).  From  Eqs. 

-1  e 

(5.14)  and  (2.7)  with  K  =  e  it  follows  that 


Under  Rayleigh  scatter  conditions,  moreover,  it  follows  from  Eq.  (2.8) 
that 


Similar  definitions  for  T||  can  be  deduced  from  Eq.  (5.15).  As  a 
practical  matter,  T||  can  be  applied  as  an  upper  bound  on  T£  when  v  ^  is 
small  as  suggested  by  L.  Wittwer.  To  accomodate  properly  geometric 
and  source  variations  in  an  extended  medium,  the  argument  of  the 
exponential  in  Eq.  (5.12)  must  be  integrated  along  the  propagation  path. 
The  integrated  equations  can  then  be  inverted  to  evaluate  or 
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VI  COHERENCE  iSANuWlDTH 


Loss  of  coherence  bandwidth  causes  group-delay  jitter  and  pulse 
dispersion.  Such  effects  are  potentially  disruptive  to  spread  spectrum 
systems  such  as  the  Global  Positioning  System  (GPS).  The  dispersive 
effects  can  be  characterized  by  the  temporal  moments 


M 


n 


f  tn  < ! vo<t> |2>  dt 


(6.1) 


Yeh  and  Liu  (1977)  have  shown,  for  example,  that 

F/2 

2nt 


T  -  M 

Td  =  M1 


F 

iS  f 


.'•F/2 


!v(f)|2[5fj  RCsfi  f  +  v]6f-odf 


(6.2) 


and 


^d  =  ^  ^  | v(t) | C  dt 


_i_  r 

2jci 

J-l 


F/2 

C  rt 

2 

|  lv(f)!2 

-2—  R(5f;  f  +  f  ) 

J-F/2 

36  f 

df 


5  f=0 


(6.3) 


where  F  is  the  bandwidth,  f£  is  the  center  frequency,  and 

R(6f;f)  =  <u(t,f  +  §f/2)u*(t;f  -  5f /2)>  (6.4) 

is  the  single- point,  two- frequency  coherence  function. 

To  calculate  R(6'f;f)  we  use  the  Huygens-Fresnel  formula  Eq.  (4.16), 
Because  Eq.  (6.4)  applies  only  to  a  single  point,  however,  it  is  un¬ 
necessary  to  preserve  the  full  generality  of  Eq.  (4.16).  We,  therefore, 
need  only  consider 
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u(p±)Re >f “  2jtiR 


(  2  ^1  )  -*  “* 

expj-i.?  2^“ju(Px  +  Sj0;^)  d? 


(6.5) 


from  which 


R(6f;f) 


1  klk2  ff  (  /£2  klk2 
2*i  (kx  -  k2)R£  JJeXP|  2  (kx  -  k2) 


z  i  Ru<^ : 

Xj 


(6.6) 


follows  after  a  variable  change  and  some  straightforward  manipulations. 
In  Eq.  (6.6)  and  are  the  wave  numbers  corresponding  to  fj  = 

f  -  6f/2  and  f^  »  f  +  6 f /2  respectively  and  pj_  is  measured  transverse 
to  the  line  of  sight. 

To  calculate  R^(Ag ;f ^jf^) j  we  use  the  phase-screen  model.  If 
e  =  fif/(2f)  is  small,  it  can  be  shown  that 


Ru0df;6f  ;f) 


exp|-D6»(^)i 


exp 


L 


9  2  2J 
2e 


(6.7) 


where  and  D(Ag)  are  evaluated  at  the  mean  frequency,  f.  (Note  that  f 
need  not  equal  f  ) . 


By  making  a  consistent  approximation  of  (k^k^) /(k^-k^) ,  changing 


variables  in  Eq.  (6.6),  and  substituting  into  Eq.  (6.7),  we  obtain  the 
general  result 


R(6  f 


-4 


exp j-i  Z£2/2.j  expj-Dg0(Z?[2Z  /f]1/2)j  dZ£  exp  2e2} 


(6.8) 


Since  we  are  working  in  a  transverse  coordinate  system,  D$0(££)  is 
functionally  dependent  on  the  quadratic  form 


C'&l  -  B^gx^y  +  A  Vg2 
A'c'  -  B/2/4 


(6.9) 
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where  A7,  B7,  and  C7  are  defined  by  Eqs.  (26a,  26b,  and  26c)  by  Rino 
(1979a).  For  convenience  the  definitions  are  summarized  in  the  Appendix. 

By  first  rotating  variables  to  remove  the  term  and  then 

performing  a  series  of  variable  changes,  Eq.  (6.8)  can  be  reduced  to  the 
single  integral 


CO 

R(Sf;f)  =  iv'7^7  J0(u»)  exp  {-iojp} 

Jo 

X  expj-D^  ^[tu2Z  6f/fl1/2)j  dtu  exp  {-a2  2e2} 

where 

A"  -  ~(A7  +  C7  +  D 7) 

C"  =  ^(A7  +  C7  -  D7) 

D7  =  v^(A7  -  C7)2  +  B7* 
a  =  (A7/  -  O'7 )  /4 
8  -  (A"  +  c")/4 


(6.10) 

(6.11a) 
(6.11b) 
(6.11c) 
(6. lid) 
(6. lie) 


and  2  is  defined  by  Eq.  (3.17).  The  same  transformation  was  used  to  derive 

Eq.  (3.14)  [cf.  Eqs,  (3.16a),  (3.16b),  and  (3.16c)],  For  isotropic 

irregularities  k"  =  c"  -  1  so  that  p  =  0.5  and  a  =  0.  For  highly  aniso- 

2 

tropic  irregularities,  k"  ~  a  where  a  is  the  axial  ratio,  and  c"  =  1. 

2  2 

For  large  a,  a  ~  p  although  p  -  a  -  k"c"/k. 

If  the  small  qQ  approximation  Eq.  (5.6)  is  used, 

D([<u2Z6f/f]1/Z)  ~  C2  [2Z6f/f]V"0,V"°'5  (6.12) 


Unfortunately  Eq.  (6.10)  cannot  be  evaluated  analytically  for  arbitrary 
v  values  within  the  admissible  range  5  <  v  <  1.5.  However,  if  v  in 

2 

F,q.  (6.12)  is  set  equal  to  1.5  without  changing  the  definition  of  Cc<i, 

o  w 

then  Eq.  (6.10)  can  be  evaluated  giving  the  result, 


(6.13) 


v/p2  “  exp  2e2} 

R(6fjf)  a  -■  - - -  rM - 

y<p  -  iHjsf/f i>2  -  a2 

where 

H  =  C2  2Z  (6.14) 

op 

Alternatively,  we  can  use  the  quadratic  approximation  Eq.  (5 #8),  whereby 

D  ([(U2Z5  f/f]1/2)  a  ^2  a2  Diq2j2z|6f/fl  (6.15) 

This  gives  rise  to  a  form  identical  to  Eq.  (6.13),  but  with  H  replaced 
by 

H,  -  2<4  Di%  <6-16> 


The  form  Eq.  (6.13)  can  be  regarded  as  exact  for  the  more  steeply  sloped 
SDFs  where  the  quadratic  approximation  can  be  applied. 


If  Eq.  (6.13)  is  used  in  Eqs.  (6.2)  and  (6.3),  assuming  that 
1 6  f /f |  «  1,  then 


(6.17) 


(6.18) 


2  2  2 

For  highly  anisotropic  irregularities  j3/(g  -q,)~1  and  0^  =  2r^. 

Equations  (6.17)  and  (6.18)  are  simplified  forms  of  similar  expressions 
derived  by  Fante  (1978)  and  Yeh  and  Liu  (1977).  The  simplifications 
arise  mainly  because  the  phase-screen  model  was  used. 
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To  give  some  idea  of  the  order  of  magnitude  of  the  expected  delay 

jitter,  in  Figure  1,  is  plotted  against  Cg.  The  parameter  values 

listed  on  the  figure  are  the  same  as  were  used  for  the  equatorial 

calculations  presented  by  Rino  (1979a).  The  quadratic  approximation  is 

also  plotted  for  i  =  q  /2jt  =  10  km.  The  results  show  that  under 
o  o 

conditions  that  produce  strong  gigahertz  scintillation,  i  ns  to  3  ns  of 
delay  jitter  can  occur,  which  is  in  agreement  with  the  calculations  of 
Yeh  and  Liu  (1977). 

The  results  have  not  yet  been  tested  against  real  data.  Analysis 
of  the  Wideband  satellite  data  is  currently  being  pursued;  moreover, 
a  GPS  receiver  system  will  be  operated  at  the  equator  in  the  summer  of 
1980,  in  part  to  evaluate  the  receiver  performance  under  disturbed 
conditions.  For  the  present,  the  model  formulas  are  adequate  for  system 
performance  evaluation.  It  should  be  noted  that  only  the  random  com¬ 
ponent  of  the  group  delay  is  included  in  Eqs.  (6.17)  and  (6.18).  The 
steady  component  must  be  added.  The  smooth  dispersive  delay  constitutes 
a  bias  that  must  be  accommodated  in  using  the  GPS  system. 
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FIGURE  1  THEORETICAL  CALCULATIONS  OF  DELAY  JITTER  CAUSED 
BY  LOSS  OF  COHERENCE  BANDWIDTH 


VII  ANGUIAR  SCATTERING 


To  analyze  radar  system  performance,  we  need  a  direct  measure  of 
the  angular  fluctuations  of  the  incident  wavefront.  As  we  have  computed 
it,  the  mutual  coherence  function  Eq.  (5.5)  gives  the  three-dimensional 
spatial  correlation  of  the  wavefield.  To  compute  the  angular  spectrum, 
therefore,  we  must  compute  the  two-dimensional  Fourier  transform  of 
Ru(Ad,Az)  in  a  plane  perpendicular  to  k. 

In  the  plane  normal  to  1c,  the  mutual  coherence  function  has  the  form 

Ru<4pA)  -  exp j-  ~D60CyX>!  (7.1) 

where 


yj.  = 


c  -  B  '^Pix^PXy  +  A 


A/C/ 


B/2/4 


1/2 


(7.2) 


The  x-axis  here  lies  in  the  plane  of  k  and  the  magnetic  field  vector. 
The  form  of  the  structure  does  not  change  aside  from  replacing  A,  B, 
and  C  in  Eq.  (3.4)  by  A7,  B7,  and  C7.  As  was  the  case  for  the  single¬ 
point  two- frequency  coherence  function,  however,  the  integral  can  only 
be  evaluated  analytically  under  the  quadratic  approximation. 


Using  Eq.  (5.10),  it  follows  from  Eq.  (7.1)  that 


^  ^(^lXJ  ^Xy) 


1 

2*DvC2 


2D  C2j 

V  > 


(7.3) 


where 

2  A'K4  +  B7KxxKXv  +  C7Kx2 
A7C7  -  B/2/4 


(7.4) 
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and  C  is  the  spherical  wavefront  correction  Eq.  (4.18).  To  convert  to 
scattering  angles,  we  need  only  replace  ICx  and  Kx  by  kfitf  and  k60 

x  y  x_,  y 

respectively.  6$,.  is  the  scattering  angle  in  the  plane  of  k  and  the 

A 

magnetic  field  vector. 

2 

Note  that  qx  “  const,  defines  a  quadratic  form.  The  axial  ratio 
is  given  by  the  formula 

AR  =*  [(A7  +  C7  +  D/)/(A/  K'  ■  D7))1/2  (7.5) 


where 

D7  ■  7 (A 7  -  C7)2  -,L  B/Z 
The  maximrm  scattering  occurs  at  an  angle 

_ _ 

A7  -  C7  . 


(7.6) 


(7.7) 


relative  to  the  x  axis. 

The  material  in  this  section  is  a  straightforward  generalization  of 
earlier  work  by  Hardin  and  Tapper t  (1974).  The  results  herein  characterize 
the  average  angular  scattering  in  a  three-dimensional  power- law  medium, 
subject  to  the  limitations  Imposed  by  the  approximations  used  to  derive 
a  simple  analytic  form  for  the  angular  spectrum. 
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Appendix 

ANISOTROPY  COEFFICIENTS 


In  Rino  and  Fremouw  (1977)  it  is  shown  that  the  two- dimensional 

phase  SDF,  $  .(ic),,  is  related  to  the  three-dimensional  irregularity  SDF 
_  60 

(K,K  )  by  the  relation 


§g^(K)  ~  r2X2  L  sec2  0  § 


£N 


.(*•  - 


tan  0 


A 

V 


(1) 


Now,  if  we  assume  that  there  exists  a  coordinate  system  (r,  s,  t)  in  which 

which  the  three-dimensional  spatial  correlation  depends  only  on 

[Ar2  +  (As/a)2  +  (Ar/r)2] with  a  and  b  2.  1,  the  corresponding 

— ♦  — > 

functional  dependence  of  on  K  and  Kz  can  be  computed  by 

straightforward  manipulations. 

The  functional  dependence  is  given  In  term s  of  the  quadratic  form 
[  (K,  Kz)TC(K,ICz)  ]1//2  where  £  is  a  3  x  3  matrix  with  elements: 


=  a2  cos2  +  sin2  ^(b2  sin2  5  +  cos2  6) 

"22  2 
C22  "  13  COS  ®  +  Sin  ® 

a  o  o  2  2  2  2 

C^3  =  a  sin  ^  +  cos  (b  sin  §  +  cos  6) 


C12  "  °21 


(d  -  1)  sin  \j(  sin  6  cos  6 


C13  “  C31  =  “  b  Sin  5  ”  COS  s*n  ^  C0S  ^ 


C23  “  C32 


(b2  -  1)  cos  ijt  sin  5  cos  6 


(2a) 

(2b) 

(2c) 

(2d) 

(2e) 

(2f) 


Note  that  if  b  =  1,  then  6  2  =  ^1  "  ^23  “  ^32  =  °*  The  an8le  t  is  the 
dip  angle.  The  angle  6  is  the  orientation  angle  of  the  transverse 
irregularity  axis.  The  angle,  0,  is  the  zenith  angle  of  k  and  c p  is  the 
magnetic  azimuth. 
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— T 


T 


The  coefficients,  A,  B,  and  C,  are  readily  deduced  by  substituting 

2 

-tan  0  atm  •  K  for  K  and  collecting  the  terms  multiplying  K!  ,  K  (C  ,  and 

^  JL  2  X  X  y 

Ky.  The  result  of  these  manipulations  gives  the  desired  expressions: 


*  2  2  A 

A  *  tan  0  cos  cp  -  2C^3  tan  0  cos  cp j  (3a) 

A  A  2 

B  =  2[C^2  +  C^g  tan  0  sin  cp  cos  cp 

A  A 

-  tan  0  (C^^  sin  9  +  C23  cos  <P)1  (3b) 

C  =  [C22  +  C32  tan2  0  sin2  cp  -  2C23  tan  0  sin  cp]  .  (3c) 


We  note  that  if  a  =  b  =  1  (isotropic  irregularities),  then 

2  2 

A  =  1  +  tan  0  cos  cp  (4a) 

2 

B  =  2  tan  0  sin  cp  cos  cp  (4b) 

2  2 

C  =  1  +  tan  0  sin  cp  .  (4c) 


Thus,  the  diffraction  pattern  is  not,  in  general,  isotropic,  even  if  the 
irregularities  themselves  are  isotropic. 

In  most  analysis  [e.g,,  Briggs  and  Parkin  (1963),  Singleton  (1973), 
Rufenach  (1975)]  a  coordinate  system  with  the  z  axis  along  the  line  of 
sight  is  used.  As  discussed  in  Rino  and  Fremouw  (1977),  this  does  not 
immediately  give  the  full  three-dimensional  structure  of  the  v/avefield. 
In  any  case,  the  appropriate  coefficients  for  the  "Briggs- Parkin"  system 
are 


A' 

B' 

c' 


2  2  2 
[A  cos  cp  +  B  cos  cp  sin  cp  +  C  sin  cp]  cos  0 

2 

[B  cos  2  cp  +  (C  -  A)  sin  cp]  cos  0 

2  cd  2 

[A  sin  cp  -  B  sin  cp  cos  “  +  C  cos  cp] 


(5a) 

(5b) 

(5c) 
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If  Eqs.  (4a),  (4b),  and  (4c)  are  substituted  into  Eqs.  (5a),  (5b), 
and  (5c),  it  follows  that  A  =  C  =  1  and  B  =  0.  Thus,  the  diffraction 
pattern  for  isotropic  irregularities  is  isotropic  in  a  plane  normal  to 
propagation  direction.  This  can,  of  course,  be  deduced  by  symmetry 
arguments . 
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